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The radiative self-energy correction to the bound-electron g factor of 2Pi/2 and 2P3/2 states in 
one-electron ions is evaluated to order a{Za)^. The contribution of high-energy virtual photons 
is treated by means of an effective Dirac equation, and the result is verified by an approach based 
on long-wavelength quantum electrodynamics. The contribution of low-energy virtual photons is 
calculated both in the velocity and in the length gauge and gauge invariance is verified explicitly. 
The results compare favorably to recently available numerical data for hydrogenlike systems with 
low nuclear charge numbers. 

PACS numbers: 12.20.Ds, 31.30.js, 31.15.-p, 06.20.Jr 



I. INTRODUCTION 

When a bound electron interacts with an external, uni- 
form and time- independent magnetic field (Zeeman ef- 
fect), the energetic degeneracy of the atomic energy levels 
with respect to the magnetic projection quantum num- 
ber is broken, and the different magnetic sublevcls split 
according to the formula 



AE = gj ijb B fj, 



(1) 



where gj is the bound-electron (Lande) g factor, ^lb = 
— e/(2TO) is the Bohr magneton, and B is the magnetic 
field which is assumed to be oriented parallel to the quan- 
tization axis. Finally, /i is the magnetic projection quan- 
tum number of the electron; i.e. the projection of its total 
angular momentum (divided by K) onto the quantization 
axis. 

In leading order, the bound-electron g factor is deter- 
mined by nonrelativistic quantum theory and is equal to 
a rational number for all bound states in a hydrogen- 
like ion. Both relativistic atomic theory as well as quan- 
tum electrodynamics (QED) predict deviations from the 
nonrelativistic result. The relativistic effects follow from 
Dirac theory and can be expressed in terms of a power 
series in the parameter Za, where Z is the nuclear charge 
number and a is the fine-structure constant. The QED 
effects are caused mainly by the anomalous magnetic mo- 
ment of the electron, which is turn in caused by the ex- 
change of high-energy virtual photons before and after 
the interaction with the external magnetic field. Here, 
by "high-energy" we refer to a virtual photon with an 
energy of the order of the electron rest mass. A second 
source for QED effects are exchanges of virtual photons 
with an energy commensurate with the atomic binding 
energy scale, which is smaller than the electron rest mass 
energy by a factor {Zaf'. Here, the electron emits and 
absorbs a virtual photon before and after the interaction 
with the external magnetic field, undergoing a virtual 
transition to a excited atomic state in the middle. For 
P states, the latter effects lead to a correction to the 
bound-electron g factor of order a{ZaY . The complete 
result for the correction of order q;(Zq;)^ is obtained af- 
ter adding the anomalous magnetic moment correction 



(high-energy part) and the low-energy photon contribu- 
tion of the same order. 

Previous studies of the bound-electron g factor for P 
states in hydrogenlike systems include Refs. [H-Ill- Quite 
recently, the problem has received renewed interest 0] . 
For few-electron ions, the bound-electron g factor has 
been investigated in Refs. [6l-[l0j. For the 2^P states of 
helium, there is still an unresolved discrepancy of theo- 
retical and experimental results (sec Refs. |10l - [l3 |V 

The expansion of the quantum electrodynamic radia- 
tive correction to the electron g factor, which is an ex- 
pansion in powers of a for a free electron, is intertwined 
with an expansion in powers of Za for a bound electron 
(this fact has been stressed in Ref. [3). For an nPj 
state in a hydrogenlike system, we can write down the 
following intertwined expansion in powers of a and Za, 



5g{nPj) = goQ + (Za) 



2 520 



0{ZaY 



+ ^Uoo + {Zaf^ + 0{ZaA. (2) 

The coefficients 1700 and 1720 characterize the relativistic 
effects, whereas 600 a-nd 620 are obtained from the one- 
loop radiative correction. The nonrelativistic result for 
the Lande g factor reads 



goo{nPi/2) 



ffoo("-^3/2) 



(3) 



The relativistic correction follows from Breit theory and 
the Dirac equation in an external magnetic field 0, , 
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g2o{nPi/2) = , 520(^^3/2) = ■ (4) 

The leading correction due to the anomalous magnetic 
reads as (see Refs. [H, Q), 



booinPi/2 



booinP3/2) = ^- (5) 



We are concerned here with the evaluation of the 620 
coefficient of nPj states, which is determined exclusively 
by self-energy type corrections (vacuum polarization does 
not contribute). 



2 



We adopt the following outline for this paper. In 
Secini we reexamine the contribution of high-energy vir- 
tual photons (see also Ref. Q). Two alternative deriva- 
tions are presented, which are based on an effective Dirac 
equation (Sec. Ill Ap and on an effective low-energy long- 
wavelength quantum electrodynamic theory (Sec. Ill B|) 
which is obtained from the fully relativistic theory by a 
combined Foldy-Wouthuysen and Power-Zienau trans- 
formation [13] • The low-energy part is also treated in 
two alternative ways. The velocity-gauge calculation in 
Sec. IIII Al ls contrasted with the length-gauge derivation 
in Sec. IIIIBI Conclusions are reserved for Sec. IIVI Nat- 
ural units (ft = c = eo = 1) are used throughout the 
paper. 



II. HIGH-ENERGY PART 
A. Effective Dirac Equation 

In Ref. |3j], the contribution to 620 due to high-energy 
virtual photons was obtained on the basis of the two- 
body Breit Hamiltonian. Here, we perform the calcula- 
tion using a simple approach, based on an effective Dirac 
Hamiltonian (see Ch. 7 of Ref. [3l ) ■ For an electron in- 
teracting with external electric and magnetic fields, this 
equation reads 



ad 



I m 



(6) 



turbative expansion of Ejj in a magnetic field read 



a ■ \ B X r 



2m 



(9) 



e K 
2m 



i,^.VV.)_ 



a - {B x r) 



where Q = \ (1—7°) is a projector onto virtual negative- 
energy states, and ip is the relativistic wave function. An 
evaluation of the first term on the right-hand side of ([9]) 
with Dirac wave functions confirms the results for 500 
and 520 given in Eqs. ([3]) and ([4]). The second term on 
the right-hand side of © yields the result for 600 as given 
in Eq. ([5]). When the Dirac wave functions are properly 
expanded in powers of Za, the second and third terms on 
the right-hand side of ([9|) yield the following high-energy 
contribution 630'^'' to the 620 coefficient defined in Eq. ^ , 



-'20 



1 
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(10) 



These results are in agreement with those given in Eq. (5) 
of Ref. [3. 



We here take into account the Dirac form factor Fi and 
the Pauli form factor i^2- The matrices a = 'j'^j and 
l3 = are the standard Dirac matrices in the Dirac rep- 
resentation (TBj . m is the electron mass, and e = — |e| is 
the electron charge. Up to the order relevant for the cur- 
rent calculation, we may approximate both form factors 
in the limit of vanishing momentum transfer as 



Fi(V2)«Fi(0) = l, i^2(V2) «F2(0) 



a 
2^ 



(7) 



The vector potential A corresponds to a uniform exter- 
nal magnetic field, i.e. A ~ ^ (^B x , and the electric 

field E is that of the Coulomb potential {eE = —W). 
Finally, V ~ —Za/r is the binding potential. So, 



-ffrad ~ a- p + 



v--d- 

2 



(^B X 



2m 



7- vy ■ 



2m 



Kpj:-B. 



(8) 



Dirac eigenstates fulfill {a-p + I3m + V)il' = -E'd'0, where 
E]j is the Dirac energy. The first few terms in the per- 



B. Long— Wavelength Quantum Electrodynamics 



It is instructive to compare the fully relativistic ap- 
proach outlined above to an effective nonrelativistic the- 
ory. In Ref. [13, a systematic procedure has been 
described in order to perform a nonrelativistic expan- 
sion of the interaction Hamiltonian for a light atomic 
system with slowly varying external electric and mag- 
netic fields. This procedure involves two steps, (i) a 
Foldy-Wouthuysen transformation of an interaction of 
the type (O, suitably generalized for many-electron sys- 
tems, and (ii) a Power-Zienau transformation to express 
the vector potentials in terms of physically observable 
field strengths. The result is an interaction, given in 
Eq. (30) of Ref. [1^, which describes a nonrelativistic 
expansion of the atom-field interaction in powers of Za 
and can be used in order to identify terms which con- 
tribute at a specified order. 

If wc are interested in evaluating the corrections to 
the g factor up to order a{Za)'^, i.e. all corrections listed 
in Eq. the relevant effective interactions for a one- 



3 



electron system are 

3 



mag 



Hi 



Hm+J2^^^ Hm = P^B{L + a)-B , (lit 



f{L + <j)-B, 



4m 



{f X a) ■ (f X B) 



= - 



2inn? 



{p-a){p-B) 



(lib) 



(11c) 



(lid) 



where /is = — e/(2m) is the Bohr magneton. We denote 
the Schrodinger-PauH two-component wave function by 
(j) in order to distinguish it from the fully relativistic wave 
function tjj. Specifically, (f) reads as (j){r) = i?(r)x^(f) 
in the coordinate representation, where R(r) is the non- 
relativistic radial component of the wave function and 
X^(f)js the standard two-component spin-angular func- 
tion |17| . An evaluation of the perturbation 



AS: 



(12) 



confirms the resuhs of Eqs. (O, g]), © and ^ for 
the high-energy part. No second-order effects need to 
be considered in this formalism up to the order in the 
Za-expansion relevant for the current study. 



III. LOW-ENERGY PART 
A. Velocity Gauge 

The most economical approach to the calculation of the 
low-energy contribution of order a{Za)'^ to the g factor of 
P states consists in a calculation of the orbital gg factor, 
with a conversion of the orbital factor to the Lande gj 
factor in a second step of the calculation. In the order 
q;(Zq;)^, one may indeed convert the spin-independent 
correction to the orbital gg factor to a spin-dependent 
correction to the gj factors of the 2Pi/2 and 2P3/2, as de- 
scribed in Appendix VK\ However, a more systematic ap- 
proach to the problem, which is also applicable to higher- 
order (in Za) corrections, is based on a perturbation of 
the nonrelativistic self-energy of the bound electron by 
the magnetic interaction Hamiltonian (jll[) . This is the 
approach outlined below. 

We thus investigate the perturbation of the nonrela- 
tivistic bound-electron self-energy |18j due to the mag- 
netic interaction Hamiltonian 



H 



M 



e 
2m 



L + (t] -B 



(13) 



given in Eq. (|lla.p . In the velocity gauge, the interaction 
of the electron with the vector potential A of the quan- 
tized electromagnetic field is given by the term —p-A/m, 



where p is the electron momentum. In an external mag- 
netic field, it is the physical momentum p— | {B xr), not 
the canonical momentum p, which couples to the quan- 
tized electromagnetic field. This amount to a correction 
SJ to the electron's transition current given by 



e 
2m 



5r = -— [B xr 



e 
2m 



e'^'^B^r^ 



(14) 



Because of the symmetry of the problem (Wigner- 
Eckhart theorem) . we may fix the axis of the B field to 
be along the quantization axis {z axis) and the projection 
of the reference state to be /i = i. This procedure al- 
lows one to simplify the angular algebra. It is inspired by 
the separation of nuclear and electronic tensors that are 
responsible for the hyperfine interaction. Such a separa- 
tion has been used in Eq. (1) of Ref. [l^ and in Eqs. (10) 
and (11) of Ref. In the case of the g factor, the 

magnetic field of the nucleus is replaced by the external, 
homogeneous magnetic field of the Zeeman effect. 

We divide out a prefactor — e/(2m) from both the mag- 
netic Hamiltonian Hm and from the current 5J^ and ob- 
tain the perturbative Hamiltonian ft.M,o and the scaled 
current ijjg. In the spherical basis, this procedure leads 
to the operators 



'A/,0 



+ o-Q : 



5jl 



(15) 



which are aligned along the quantization axis of the ex- 
ternal magnetic field. The index zero of the operators, 
in the spherical basis, denotes the z component in the 
Cartesian basis (see Ref. [13 )■ The following shorthand 
notation for the atomic states with magnetic projection 
= i proves useful, 

EE \nP,{^l = i)) , {hMfl) EE {]\\hMM) ■ (16) 

Finally, we can proceed to the calculation of the per- 
turbed self-energy. The nonrelativistic (Schrodinger) 
Hamiltonian of the atom is 



NR. 



2m 



(17) 



and the nonrelativistic self-energy reads 



SE= - 



2a 



diouj 



3h 



H^K — E] 



NR. 



(18) 



The wave-function correction to the self-energy reads 
4a 



SE. 



dujuj 



(19) 



1 



^^NR — ^'nR, + ^ 



1 



E^vt — H 



NR. 



Here, [1/(£'nr — ^^nr)]' is the reduced Green function, 
with the reference state being excluded from the sum over 
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2a 



1 



In 



(Za) 



0.25134j , 
(25a) 



0.15907 

(25b) 



virtual states. The contribution of virtual nPj states A numerical evaluation leads to the following results, 
(with j being equal to that of the reference state and 
n > 2) vanishes because of the orthogonality of the non- 
relativistic radial wave functions. The interaction /ij\/,o 
couples 2Pi/2 and 2P3/2 states, but the contribution of 
virtual states with different j as compared to the refer- 
ence state vanishes after angular algebra [l3| because the 
self-energy interaction operator {p[l/(_ffNR~£'NR+w)]p} 
is diagonal in the total angular momentum. Virtual 
states with different orbital angular momentum than the 
reference state are not coupled at all to the reference state 
by the action of the perturbative Hamiltonian /im,o- Be- 
cause all contributions vanish individually, we can thus 
conclude that SEfj, = 0. 

Hence, we have to evaluate first-order corrections to 
the Hamiltonian and to the energy corresponding to the 
replacements -ffNR, Hnr + ^a/,o and Snr ^ -Enr + 
(^j\/,o) in Eq. (jl8p . Furthermore, we have a correction 
to the current corresponding to p/m -> p/m + 6jo- The 
energy correction reads as 



The correction to the current is given by 



SEc = - 



X U 



1 



P 



Hnr — -Enr + w 



(26) 



SEe 



2a 



duj L 



where we take into account the multiplicity factor due to 
the current acting on both sides of the propagator. The 
corresponding correction to the g factor is 



{hM,o) 



(20) 



6gc{nPj) ^ gooinPj) 



SE, 



Ht^r — E 



NR 



where me is an upper cutoff for the photon energy [l8| 
(the scale-separation parameter e is dimensionless) . It 
corresponds to the following g factor correction. 



We obtain the following numerical results, 

2a 

Sgc{2Pi/2) = :^(^a)' 0.24607, 



SgsiriPj) = gm{nPj) 



6Ee 
{hM,o) 



5gc{2Py2) 



(21) 



4q; 
3^ 



{Zaf 0.06151, 



(27) 



(28a) 



(28b) 



A numerical evaluation of this correction according to 
established techniques [2l| yields 



„(2P3,,)^|(Za)^(-il„(^ 



Summing all low-energy corrections, the spurious loga- 
rithmic terms cancel, and we obtain 



4o'(2Pi/2)= 0.98428, 



0.12831 j , 
(22a) 



0.12831 ] 



4o^(2i^3/2) 



0.49216, 



(29a) 
(29b) 



(22b) 



For the correction to the Hamiltonian, we get 
2a 



SE 



H 



377771^ Jq 



1 



Hnr — Ej<sji + Lo 



h-MflP 



(23) 



where we have used the relation [Hnr, ^j\f,o] = 0. This 
translates into the following correction for the g factor. 



as the spin-dependent low-energy contribution to the 
bound-electron g factor. The result for 2P3/2 is equal to 
half the correction for 2Pi/2 (this fact is independently 
proven also Appendix|X|. We denote the low-energy con- 
tribution to the 620 coefficient defined in Eq. ^ by 620'' ■ 



B. Length Gauge 

In the length gauge, the interaction with the quan- 
tized electromagnetic field is given by the dipole interac- 
tion —ex-E, where E is the electric-field operator p^ . 
The gauge-invariant [22| nonrelativistic self-energy in the 
length-gauge reads 



6E 



SgniriPj) = goo{nPj) 



SE 



H 



(24) 



2a 



dlOiO^ 



3\ 



1 



^^NR — E} 



NR, 



(30) 
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In the length gauge, the contribution of the wave-function 
correction vanishes because of the same reasons as for the 
velocity gauge. Also, there is no correction to the transi- 
tion current, because the canonical momentum does not 
enter the interaction Hamiltonian in the length gauge. 
We only have corrections to the Hamiltonian and to the 
energy. We start with the energy perturbation. 



6-K 



duj uj 



(31) 



Hi 



NR - 



NR 



The subscript £ instead of E serves to differentiate the 
length-gauge as opposed to the velocity-gauge form of 
the correction. Indeed, the numerical results for the cor- 
responding correction to the bound-electron g factor are 
different from those given in Eq. ([22) and read 



- 0.88488 ) , 



,,,(2P3,,,.g(Z„,= (4,„(_ip 



For the correction to the Hamiltonian, we get 
6 En ^ — I dujuj^ 



(32a) 



(32b) 



Stt Jo 



1 



Hnr — -Bnr + w 



U/,0 X 



A numerical evaluations leads to 
S9nm/2) = ^{Zaf In 



3tt' ' \2"'\ {Zaf 



(33) 



1.25399j , 

(34a) 



2-V(Z.)V°™j- 
(34b) 



Adding the length-gauge corrections, the logarithmic 
terms cancel, and it is straightforward to numerically ver- 
ify the gauge-invariance relation 



and with the help of a somewhat lengthy calculation, it 
is possible to show analytically that the velocity- gauge 
and the length-gauge forms of the low-energy contribu- 
tions are equal. The calculation follows ideas outlined in 
detail in Ref. [l^ where the more complicated case of a 
relativistic correction to a transition matrix element was 
considered. Here, we are interested mainly in the numeri- 
cal value of the correction, for which the gauge invariance 
provides a highly nontrivial check. Note that the matrix 
elements governing the transitions of the reference to the 
virtual states are completely different in the length and 
in the velocity gauges, and the final results are obtained 
after summing over the discrete and continuous parts of 
the spectrum of virtual states. 



IV. CONCLUSIONS 

In our approach to the calculation of the bound- 
electron g factor of P states, the contribution due to 
high-energy virtual photons can be obtained using two 
alternative approaches, based either on an effective Dirac 
equation or on a low-energy effective Hamiltonian. The 
contribution due to low-energy photons is treated as a 
perturbation of the bound-electron self-energy [l^ due 
to the interaction with the external uniform magnetic 
field. Corrections to the Hamiltonian, to the bound-state 
energy and (in the velocity gauge) to the transition cur- 
rent have to be considered. The final results for the low- 
energy parts in the velocity- and length-gauges agree al- 
though the individual contributions differ (including the 
coefficients of spurious logarithmic terms). 

Adding the high-energy contribution to the g factor 
correction given in Eq. (|10p and the low-energy effect 
given in Eq. ((29)) . we obtain the following results for 
the self-energy correction of order a{Za)'^ to the bound- 



electron g factor of 2P states (620 = ^20^'' + ^20 
fe2o(2Pi/2) = 0.48429, 



(37a) 



&2o(2P3/2) = 0.59214, (37b) 

where the 620 coefficient has been defined in Eq. Both 
above results compare favorably with recently obtained 
numerical data for low-Z hydrogenlike ions [23] (see also 
Appendix|B|) . An obvious generalization of the formalism 
outlined here to the 3P and AP states yields the results 

620(3^1/2) = 0.40500, 



^20(3^3/2) = 0.55250, 



(37c) 
(37d) 



5gE + 5gH + 5gc = Sgs + Sg-u 



(35) 



&2o(4Pi/2) = 0.31331, 



(37e) 



and thus, the numerical results already given in Eq. (|29p . 

Let us finally discuss the analytic proof of the gauge 
invariance. Using the commutator relation 



PL 
m 



(36) 



620(4^3/2) = 0.50665. 



(37f) 



The two main results of the current investigation can 
be summarized as follows. First, in Sec. IIIII we formu- 
late a generalizable procedure for the calculation of low- 
energy corrections to the Lande gj factors in one-electron 
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ions, applicable to P states and states with higher angu- 
lar momenta. This procedure is based on choosing a spe- 
cific reference axis for the external magnetic field. In the 
future, it might be applied to include higher-order terms 
from the Hamiltonian ([TT|) which couple the orbital and 
spin degrees of freedom. Second, we resolve the discrep- 
ancy reported in Ref. [23| regarding the low-Z limit of 
the a (Za)^ correction to the gj factor with previous re- 
sults reported in Ref. [1| for this correction [see Eq. ((37)) 
and Appendix |B] . In Appendix \^ it is shown that the 
discrepancy to the results of Ref. ^ can be traced to 
the final evaluation of the logarithmic sums over virtual 
states, while the angular momentum algebra is in agree- 
ment. That is a further reason why the cross-check of our 
calculation in the length and velocity gauges appeared to 
be useful. 

Regarding the experimental usefulness of the obtained 
results, we can say that recent proposals [2^ concerning 
measurements of the bound-electron g factor of low-Z 
hydrogenlike ions are based on double-resonance schemes 
that also involve transitions to 2P states in the presence 
of the strong magnetic fields of Penning traps. In order 
to fine-tune the double-resonance setup, the results ob- 
tained here might be useful. Also, the results reported 
here serve as a general verification for the analytic formal- 
ism used in the theoretical treatment of corrections to 
the g factor of states in helium, for which an inter- 
esting discrepancy of experimental and theoretical results 
persists (see Ref. [13 and Sec. V of Ref. [3). 

We conclude with the following remark. Our calcu- 
lation concerns the g factor of P states, and we con- 
firm that in the order a{Za)'^, low-energy virtual pho- 
tons yield an important contribution. From the treat- 
ment in Sec. IIIIl we can understand physically why 
there is no such low-energy effect of order a{Za)'^ for S 
states. Namely, a Schrodinger-Pauli S1/2 state happens 
to be an eigenstate of the Hamiltonian /im,o- We have 
hM,o\nSl/2i^Ji = ±i)) = ±|nS'i/2(/i = ±^)). This prop- 
erty holds because an S state carries no orbital angular 
momentum, and therefore is an eigenstate of the third 
component of the spin operator ctq, and 620 therefore 
vanishes for S states [1^. For P states and states with 
higher orbital angular momenta, the situation is different: 
these states arc not eigenstate of h m,o irrespective of their 
angular momentum projection, even though ft.j\/,o com- 
mutes with the nonrelativistic Hamiltonian -/Jnr- There- 
fore, there is a residual effect of order a{Za)'^ due to 
low-energy virtual photons for states with nonvanishing 
angular momenta. 
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TABLE I: Higher-order remainder functions 41/2(^0;) and 
13/2 (-^c) for the self-energy correction to the g factor of 
2P1/2 and 2P3/2 states, respectively, as obtained recently 
in Ref. [23]. The value of a employed in the calculation is 
— 137.036, and the numerical uncertainty of the all-order 
(in Za) calculation due to the finite number of integration 
points in the numerical calculation is indicated in brackets. 



z 


9i/2{Za.) 


93/2{Za) 


1 


0.121258 (21) 


0.148104 (21) 


2 


0.121715 (22) 


0.148294 (24) 


3 


0.122414(19) 


0.148567(24) 


4 


0.123280(14) 


0.148851 (22) 


5 


0.124305(10) 


0.149338(18) 


6 


0.125473(7) 


0.149816(14) 


7 


0.126803 (5) 


0.150350(11) 


8 


0.128186 (4) 


0.150933 (7) 


9 


0.129711 (2) 


0.151561 (4) 


10 


0.131336(2) 


0.152231 (4) 
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FIG. 1: (Color online) The higher-order remainder function 
ii/2{Za) is shown as a function of Z. Numerical values for 
ii/2{Za) are given in Tabled The point at Z = is given by 
the coefficient i62o(2Pi/2) and is approached smoothly. 




FIG. 2: (Color online) Same as Fig. [T] but for the 2P3/2 state. 
The higher-order remainder function 13/2(^0=) is plotted as a 
function of Z, with numerical values for 13/2 (-^q) given in Ta- 
blelH The point at Z = is limza^o 13/2 (^q) ~ jb2a(2P:j/2)- 
The limit is approached smoothly. 



7 



Appendix A: Remarks on the Low— Energy Part 



First of all, let us remark that our results reported 
in Sec. IIIII can be expressed as logarithmic sums over 
the spectrum of atomic hydrogen. After performing the 
angular algebra [l7|. we can write the total low-energy 
correction 5g^'^\nPj) = Sg£{nPj) + 5g-u{nPj) as follows, 



X In I I - — 2_^iE,nD - Enp) 



X {nP\\r\\nDf In 



2 / 2|i?m£) — Enp\ 



(Al) 



and 



^3/2) = P E(^™5 - E^pf {mSMnPf 



97r 



X In i ^\ErnS^E^p\ ] _ 4a ^^^^^^^ _ 



{ZaYm 



97r 



/ mm n\2 1 ('^\EmD—Enp\ 

X {nP\\r\\nD) In 



(A2) 



Here, (TOS'||7^|nP) and (m£)| |r| | nP) are reduced matrix 
elements in the notation of Ref. [23|- Because {EmS — 
Enp) oc [ZaY and (m5||r||?iP) oc {Za)~^, the above 
corrections to the g factor are manifestly of order a{Zay' . 
The sums over m extend over both the discrete as well 
as the continuous part of the hydrogen spectrum and can 
conveniently be evaluated using basis-set techniques [1^ . 

Because it may not be completely evident from the 
presentation in Ref. [Ij, we reemphasize here that the 
authors of the cited article evaluate a correction 5gi to 
the orbital g factor according to the definition 



-ffoft = gs fJ-B S ■ B + gi^iB L ■ B 



(A3) 



for the effective interaction of a bound electron with an 
external magnetic field {S = measures the electron 
spin, and we have gs « 2, w 1). Here, an inter- 
esting analogy to the description hyperfine splitting can 
be drawn, because the interaction with the nuclear mag- 
netic field can also be separated into distinct components, 
namely the orbital, spin-dipole, and Fermi contact terms 
given, e.g., in Eqs. (24)— (33) of Ref. H. 

Returning to the discussion of the g factor, we see that 
as a combination of g^ and gc, the Lande g factor is ob- 
tained as 



9j = ge ■ 



j(j + i) + e{e + i)-.s{s + i) 

+ 1) 

jij + l) + s{s + l)-£i£+l) 



3j{j + i)~e{e + i) + s{s + i) 

2j(j + 1) 



(A4) 



In our notation, the analytic result given in Eq. (14) of 
Ref. for the correction to the orbital g factor reads 



SgiinP) - — E(S,nS - E^pf {mS\\r\\nPf 

■ ''^{EmD — EnpY 



X In 



2|-Em5 ~ Enp\ 



{Za)^m 



a 
3^ 



X {nP\\r\\nD)^ In 



2\E, 



mD 



En} 



{ZaYm 



(A5) 



H) read | for 



The prefactors multiplying g^ in Eq 
Pi/2 and I for P3/2 states, and therefore the analytic 
formula obtained in Eq. (14) of Ref. Q is in agreement 
with our approach for both P1/2 a-nd P3/2 states, after 
the correction to gi is converted into the correspond- 
ing modification of gj. However, their numerical result 
6gi — — 0.24q;'^ disagrees with our result both in sign and 
in magnitude. Indeed, expressed in terms of our 620'' co- 
efficient, the results indicated in Ref. Q would imply that 



40^(2/^1/2) 



-4.02 and 4o^(2P3/2) 



-8.04. 



Finally, let us note that the calculation of the orbital 
correction to the gj factor is only applicable to order 
Q!(Zq!)^, not a{Za)'^^ because the higher-order terms in 
the magnetic interaction (jlip couple the orbital and spin 
degrees of freedom. The formalism outlined in Sec. IIIII 
generalizes easily to the calculation of higher-order cor- 
rections that couple spin and orbital angular momentum, 
which might be needed in the future, whereas the sepa- 
ration into spin and orbital g factors only holds up to 
order a{Zaf'. 



Appendix B: Comparison to Numerical Data 

We parameterize the one-loop self-energy correction 
iJ^^^g to the g factor of P states as 

<5(i)5(2P,) = ^ {600 + [Zaf i,[Za)) , (Bl) 

where ij{Za) is the nonperturbative (in Za) remain- 
der function. The remainder functions ii/2{Za) and 
^3/2(^0;) for 2Pi/2 and 2P3 /^, respectively, have recently 



been evaluated in Ref. [2J]. Numerical values of ij{Za) 
for Z = 1, . . . , 10 are given in Table [J Note that these 
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data imply a spin-dependence of the higher-order correc- 
tion term ^20 beyond the spin-dependence of the high- 
energy part given by Eq. (|T0|) . The hmit as Za — > of 
the remainder function ij{Za) is 

\un^^,iZa) ^ lb,„{2P,) . (B2) 



As shown in Figs. [T] and [2l this hmit is being consistently 
approached by the numerical data. 
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